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Abstract 

We consider a system of coupled Gross-Pitaevskii (GP) equations describing a bi¬ 
nary quasi-one-dimensional Bose-Einstein condensate (BEG) with intrinsic time- 
dependent attractive interactions, placed in a time-dependent expulsive parabolic 
potential, in a special case when the system is integrable (a deformed Manakov’s 
system). Since the nonlinearity in the integrable system which represents binary 
attractive interactions exponentially decays with time, solitons are also subject to 
decay. Nevertheless, it is shown that the robustness of bright solitons can be en¬ 
hanced in this system, making their respective lifetime longer, by matching the time 
dependence of the interaction strength (adjusted with the help of the Feshbach- 
resonance management) to the time modulation of the strength of the parabolic 
potential. The analytical results, and their stability, are corroborated by numerical 
simulations. In particular, we demonstrate that the addition of random noise does 
not impact the stability of the solitons. 
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1 Introduction 


The experimental realization of Bose-Einstein condensates (BECs) in dilute 
gases [1] has provided a deep insight into the realm of macroscopic quantum 
phenomena, and has now become an experimental testbed for direct manip¬ 
ulations of matter waves. The creation of bright [ 2 ] and dark | 3 ] solitons in 
BEC has drawn additional strong interest to this area |1]. The behavior of 
the single-component (scalar) BEC is affected by the external trapping po¬ 
tential and the interatomic collisions, characterized by the scattering length. 
Actually, the lifetime of the self-attractive scalar BEC is very short, as it is 
subject to decay, once its density exceeds a critical value. Therefore, effective 
stabilization of self-attractive condensates remains an important problem. In 
this context, the use of binary (two-component) BECs, composed of either 
hyperhne states of the same atom 06], or different atomic species m , may 
be very relevant, as the dynamics of binary condensates can be controlled by 
adjusting both the intra- and inter-species scattering lengths via the Feshbach 
resonance (FR), see book |H| and references therein. In particular, the concept 
of the exchange of energy between the two components of a binary system, 
which is known in the realm of optical solitons [9], has been demonstrated in 
binary condensates as well [10]. This mechanism can be employed to increase 
the lifetime of the quasi-one-dimensional BEC. 

Recently, it was shown m that one can manipulate trajectories of bright 
solitons and their intensity distribution in the coupled nonlinear Schrodinger 
(NLS) equations by adjusting the respective self-phase-modulation (SPM) and 
cross-phase modulation (XPM) interactions. Experiments, starting from early 
works [S] , and including more recent ones 0 , made use of the repulsive inter¬ 
actions in multi-component BECs composed of ®^Rb atoms to demonstrate the 
possibility of producing long-lived multiple-condensate states. Those results, 
and the availability of multi-component condensates with attractive interac¬ 
tions [T21I13II14] , suggest one to look for a mechanism for enhancing the stability 
of multi-component BECs by engineering properties of bright solitons, which 
are coherent condensates by themselves. 

In this paper, we consider the dynamics of a binary condensate in an ex¬ 
pulsive time-dependent parabolic potential, governed by a system of coupled 
Gross-Pitaevskii (GP) equations. The motivation to study the properties of 
the condensates in expulsive potentials stems from the fact that BEGs are 
more stable in the conhning trap while they get compressed and are more 
unstable in the expulsive trap. Expulsive parabolic potentials are also used 
in experiments for probing various dynamical properties of the condensates 
[TB] . We then study the evolution of the condensates assuming that the scat- 
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tering length can be made time-dependent throngh the FR, so as to cast 
the system into an integrable form. We thus observe that the condensates 
in the expulsive temporarily transient harmonic potential may sustain stabil¬ 
ity, while their counterparts in the time-independent potential would rapidly 
decay. The analytical results are conhrmed by numerical simulations employ¬ 
ing the split-step Crank-Nicolson method. Furthermore, we demonstrate that 
the addition of noise does not impact the stability of the condensates in the 
expulsive transient potential. 


2 The model and the Lax pair 


Considering a two-component BEC, with equal atomic masses m and attrac¬ 
tive interactions (such as the condensate composed of two hyperhne states of 
^Li da or 85Rb [miT3] atoms), trapped in a parabolic potential, the mean- 
held evolution of the setting is governed by coupled quasi-one-dimensional GP 
equations im, written in a scaled form: 


.d'lpi _ 
* dt ~ 



(~2^ + + &i2|V^2p + V’l, 

(“ 2 ^ + + &22|^2P + y^^^ ^2, 


( 1 ) 


where y {j = 1 , 2 ) is the macroscopic wave function of the j-th component 
subject to the normalization condition = 1 and |V’ 2 pdx = 

N 2 /N 1 . The interaction between the atoms is described by the self-interaction 
coefficients, bjj = 4njjNi/l±, and the interaction between different compo¬ 
nents is controlled by bjk = 4:njkNi/l±, where rijj is the scattering lengths of 
the j-th components, and rijk is the same for collisions between j and k . The 
dynamics of the two-component BECs composed of hyperhne states of ®^Rb 
atoms has been experimentally investigated in detail, and it it was shown how 
the binary condensates generate various patterns m- The system of coupled 
GP equations for the two-component BEG with unequal scattering lengths is 
nonintegrable, and has been investigated earlier [IH]- In the present work, we 
consider the condensates with symmetric interaction strengths, [I01I22] assum¬ 
ing bii = 621 = bi 2 = 622 (fhe scattering lengths in attractive binary mixtures 
may be indeed made nearly equal [T^19lll5j ). and A|(t) = a;|(t)/a;^, where ujj 
and a;_L represent the angular frequencies of the trapping potential in the axial 
and radial directions, respectively. Time t and coordinate x are measured in 
units of 2 /a;j_ and /j_ = ^Jh/mu±, respectively, while /j_ represents the lin¬ 
ear oscillator length of the tight trapping potential acting in the transverse 
direction. 
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Fig. 1. Upper Panels (a-b): The analytical density profile of the bright solitons 
in the time-independent expulsive trap for F = —0.25 and g{t) = Lower 

Panels (c-d): The numerically simulated density profile for the same case. 

Further, assuming Ui{t) = U 2 it) = u{t) (i.e., Ai = A 2 = A), and allowing the 
scattering lengths bj^ and the strength of trapping potential A^ to vary with 
time, eq. ([^ takes the following form; 


+ '^ixx + 2g{t){\^jJl\'^ - 1 - = 0 

#2t + '02XX + 2g{t){\^jJl\^ + \^|J2\‘^)^|J2 + X^{t)x^il)2 = 0. (2) 


It is relevant to stress, in passing, that we consider the GP equations in the 
framework of the mean-held approximation, thus completely neglecting quan¬ 
tum correlations which, beyond the limits of the mean-held theory, may cause 
macroscopic entanglement between matter-wave solitons via their collisions 



In eqs. ([^, the nonlinearity takes the Manakov’s form [ 23 ] (which is a nec¬ 
essary condition if an integrable system is sought for), with g{t) standing for 
the common interaction strength, while \{t) represents the time-dependent 
trap frequency. Below, we consider values g{t) > 0, which correspond to the 
attractive signs of intraspecies and interspecies interactions. Under the special 
integrability condition imposed on g{t) and A(f), (see eq. ([^ below), eqs. ([^ 
admit a representation in the form of an eigenvalue problem. 
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+ = 0 , 

+ V"$ = 0, 

where $ = (01,02,03)"^ and 

Kit) Qi Q2 
-Ql -tC{t) 0 

-g^ 0 -<(t) 

nil ni 2 ni 3 ^ 

V21 V22 V23 
^31 V32 V33 

with 

nil = + *r(t)a;C(t) + -QiQl + -Q2Q2 

V12 = T{t)xQi - C{t)Qi + ^Qix 

ni3 = r(t)a;g2 — C{t)Q2 + 2 ^ 2 ® 
n 2 i = —r(t)a;g^ + C{i)Qi + 2^1* 

V22 = iCity - iT{t)xCit) - ^QiQl 
n23 = —^Q2Qi 

v,, = -T{t)xQ; + at)Q;+'-Ql, 

V32 = —-^QiQ*2 

V33 = iC{tf - iT(t)xC{t) - ^Q2Q*2 

gi = ^^0i(a;,t)e*rW"'/2 
g2 = ^^02(a:,t)e'^W"'/2 




( 3 ) 


(4) 


(5) 
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The compatibility condition, ($a;)i = leads to the zero curvature equa¬ 

tion Ut — Vx + \U,V] = 0 which yields the integrable system of coupled GP 
equations (j^, provided that the spectral parameter C,{t) obeys the following 
nonisospectral condition: 


c(t) = /xexp 


T(t)dt) , 


( 6 ) 


where /i is a hidden complex constant and r(t) is an arbitrary function of 
time, which is related to the trapping frequency by the following constraint: 

X\t) = T^{t) -dT{t)/dt. (7) 

Then, it follows from eq. Q that frequency X{t) is related to scattering length 
g{t) through the integrability condition: 


- 9 {t)d^g{t)/dt‘^ + 2 [dg{t)/dtf - X‘^{t)g^{t) = 0. (8) 

Thus, the coupled GP equations (|^ is completely integrable if the trapping 
frequency, A(t), and scattering length, g{t), are subject to constraint ([^. For 
the time-independent trap, X{t) = const = Ci, eq. ([^ yields g{t) = [2^ . 

In this work, we focus on the consideration of the integrable system satisfying 
condition governed by eq. ([^. If it is slightly broken, the result depends on the 
accumulation of the deviation from the integrability over the time interval, T, 
corresponding to essential dynamical regimes, which may be easily identihed in 
all examples displayed below. Namely, if the deviation from the integrability 
is characterized by difference AX{t) from the value imposed by eq. (|^, the 
condition for the system to remain close to the integrability is 

[ AX{t)dt -C 1. (9) 

Jo 

It should be mentioned that one can convert the coupled GP equation into 
the celebrated Manakov’s model by a suitable transformation, as shown, e.g., 
in ref. [10]. Nevertheless, the direct formulation of the integrability formalism 
in the form of eqs. 0 . 0 . and ([^ is quite useful, as it makes it possible to 
apply the gauge-transformation method for generating multisoliton solutions, 
and it may also be used for the search of more general integrable systems. 


3 Analytical and numerical results for two-component bright soli- 
tons in the integrable system 


Using the gauge-transformation approach [23], bright solitons of the coupled 
GP equation ([^, subject to the integrability condition ([^, can be found in 
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Fig. 2. Upper panel (a-b): Analytical results for the compression of the condensate 
in the time-independent expulsive trap for F = —0.5 and g{t) = 0.5e“^'^*. Lower 
Panel (c-d): Numerically simulated density profile for the same case. 

the form of 


V'S‘' = (10) 

where 


6i = 2f3ix + 4 y aiPidt — 26i, 

= 2aix + 2 y(«i - f^i)dt - 2xi, 

with «! = aio exp [/r(f)(if], /3i = /3io exp [/r(f)(if] while and xi cire ar¬ 
bitrary parameters, and are conpling coefficients, which are snbject 

to constraint + 1^2^^ P = 1- The gange-transformation approach can be 

extended to generate mnlti-soliton solntions [22] • 

Thns, it is obvions from the above that the amplitnde of the bright solitons 
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Fig. 3. Upper panels (a-b): Onset of the decay (analytical results) for F = —0.9 
and g{t) = 0.5e“^'®* in the expulsive time-independent trap. Lower Panels (c-d): 
Numerically simulated density prohle for the same case, showing the decay. 

is determined by the temporally modulated scattering length g(t) and trap 
frequency T(t) {/3i(t) varies exponentially with r(t)). 

To start the analysis of particular solutions relevant for the physical realization 
of the system, we now switch off the time dependence of the harmonic trap 
and adopt the scattering length in the form of g{t) = for which 

eqs. ([^ and ([^ render the parabolic potential expulsive, T = —0.25. The 
corresponding density prohle of the condensate is shown in the upper panel of 
fig- 0 (figs- [3a) and[3b)). The corresponding density prohle of the condensates 
produced by numerical simulations of the real time propagation, using the 
split-step Crank-Nicolson method is shown in hgs. [^c) and 0 (d). From hgs. 
0:a)0d), we observe perfect agreement between the analytical and numerical 
results. When we double the trapping frequency to T = —0.5, keeping the 
trap expulsive and choose the scattering length as g{t) = 0.5e“°'®* consistent 
with eqs.(0) and (j^, the compression of the condensates sets in, as shown 
in hgs. [^a-b). This is conhrmed by numerical simulations shown in hgs. j^c) 
andj^d). When we further enhance the trap strength, to 3.6 times the original 
value, T = —0.9, keeping the trap expulsive and choose the scattering length 
as g{t) = 0.5e“°'®* (again consistent with the integrability condition given by 
eqs. (j^ and (j^), one observes decay (spreading out) of the condensates, as 
shown in hgs. Again, results of numerical simulations, shown in hgs. 
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Fig. 4. Upper panel (a-b): The analytical density profile produced by switching on 
the time dependence of the expulsive trap for r(t) = —0.25f and g{t) = . 

Lower Panel (c-d): Numerically simulated density profile for T(t) = —0.25t and 
g{t) = 

[^c) and I (d), match with figs. ia) and|^b), respectively. 


To enhance the stability of the condensates, we now switch ON the time 
dependence of the harmonic trap, keeping it expnlsive, as T{t) = —0.25t, 
and choose g{t) = which is consistent with eqs. Q and The 

corresponding density profile is shown in figs. |^a,b). This analytical solntion 
is confirmed by the nnmerical simulations, as shown in figs. |^c,d). When we 
increase the expulsive-potential strength by a factor of 20 (as compared to 
fig.(|4)), accordingly setting g{t) = which is consistent with eqs. ([^ 

and (^, the corresponding density profile shown in figsj^(a,b) again matches 
with its numerical counterpart shown in figs, [^c) andj^d). 


When we further increase the-time dependent expulsive-trap frequency r(f) 
by a factor of 100 (as compared to fig.(|^), and choose the interaction strength 
in accordance with eqs.(I^ and ([^), an abrupt increase in the density is no 
more observed, as figs. |Ma-b) show. In fact, any further increase of T{t), 
and the respective change of the interaction strength g{t) consistent with 
eqs. (j^ and (|^, does not lead to a significant increase in the density, even 
though the attractive-interaction strength increases rapidly. In other words, 
the condensates with the attractive interactions in the time-dependent ex- 
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Fig. 5. Upper panels (a-b): Analytically found density profiles of the condensates 
in the time-dependent expulsive potential for T{t) = —5t and g{t) = 0.5e“^'®* . 
Lower panels (c-d): Numerically found density profiles for the same case. 

pulsive parabolic potential remain essentially stable over a reasonably large 
interval of time. Again,analytical results match the numerical simulations, as 
shown in hgs. g(c-d). 

Thus, we observe that the two-component BEC with attractive interactions in 
the time-dependent expulsive trap, described by the integrable system, is more 
long-lived, as its lifetime can be increased by means of the FR management, 
in comparison with the condensates in the time-independent expulsive trap. 
Note that the scalar (single-component) condensate, stabilized by means of 
the Feshbach management in a similar setting, quickly decays in the course of 
the evolution. 

To further conhrm the stabilization, we added random noise to simulations 
performed in the time-dependent expulsive trap. The respective numerically 
simulated density prohles demonstrate, in hgs. 0-0 , that the noise does 
not impact the stability of the condensates. In all the above numerical simula¬ 
tions, the norm of the wavefunctions is conserved. This conclusion corroborates 
that the two-component condensate in the time-dependent expulsive harmonic 
trap, governed by the specially devised integrable system, is more stable in 
comparison with its counterpart in the time-independent trap. The above re¬ 
sults indicate the possibility of increasing the life span of the two-component 
BEC with attractive interactions in the time-dependent parabolic potential 
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Fig. 6. Upper panels (a-b): The analytically found density profile in the time-de- 
pendent expulsive potential for r(t) = —25t and g{t) = . Lower panels 

(c-d): The numerically simulated density profile for the same case. 



Fig. 7. The numerically simulated density profile, corresponding to fig.|^ but withthe 
addition of random noise. 

experimentally, employing the FR management, which may be applied to the 
condensates composed of [13], [T3| and ^Li [231 atoms. 


4 Conclusion 


In this paper, we have shown that one can make use of the precise control of 
the scattering length by means of the FR (Feshbach resonance) to design con- 
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Fig. 8. The numerically simulated density profile, corresponding to fig. (§ , but with 
the addition of random noise. 



Fig. 9. The numerically simulated density profile, corresponding to hgj^ but with 
the addition of random noise. 

ditions under which the mean-held dynamics of two-component condensates 
in the time-dependent expulsive parabolic potential is governed by the inte- 
grable model (a deformation of the Manakov’s system, with a time-dependent 
scattering length). A regime may be selected, in which solitons stay effectively 
stable against the decay for a reasonably large interval of time, compared 
to the condensate in the time-independent expulsive parabolic potential. The 
analytical results, and their stability, have been corroborated by the com¬ 
parison with numerical simulations. In particular, it was demonstrated that 
binary condensate with the attractive interactions remains stable in the time- 
dependent harmonic potential in spite of the addition of random noise. Thus, 
the conclusion is that the bright solitons, which would quickly decay (possibly 
following short-time compression) otherwise, may be made more robust by the 
proper selection of scattering length through FR management in combination 
with the appropriate choice of the time dependence of the parabolic potential. 
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